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ABSTRACT 

The first generation of stars (commonly known as population III) are expected to form 
in low-mass protogalaxies in which molecular hydrogen is the dominant coolant. Ra- 
diation from these stars will rapidly build up an extragalactic ultraviolet background 
capable of photodissociating H2, and it is widely believed that this background will 
suppress further star formation in low-mass systems. 

However, star formation will also produce an extragalactic X-ray background. 
This X-ray background, by increasing the fractional ionization of protogalactic gas, 
promotes H2 formation and reduces the effectiveness of ultraviolet feedback. 

In this paper, we examine which of these backgrounds has the dominant effect. 
Using a simple model for the growth of the UV and X-ray backgrounds, together with 
a detailed one-dimensional model of protogalactic chemical evolution, we examine 
the effects of the X-ray backgrounds produced by a number of likely source models. 
We show that in several cases, the resulting X-ray background is strong enough to 
offset UV photodissociation in large H2-coolcd protogalaxies. On the other hand, small 
protogalaxies (those with virial temperatures T v - a < 2000 K) remain dominated by the 
UV background in all of the models we examine. 

We also briefly investigate the effects of the X-ray background upon the thermal 
and chemical evolution of the diffuse IGM. 
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1 INTRODUCTION 

In cosmological models based on cold dark matter (CDM), 
the first stars are believed to form within sm all protogalax- 
ie s, with virial temperatures T v p < 10 4 K ( Couchman & 



Rees 



1986 



Haiman et al. 



L996a; Pegmark et al. 1997). Cool 



ing within these protogalaxies is dominated by molecular 
hydrogen, IU, which fumis via Llie gas-phase leacliuns 



H + e 
H" +H 

and 

H + H+ 
H+ + H 



H~ + 7 
H 2 +c, 



H 2 +H 4 



(1) 

(2) 



(•3) 
(4) 



even in the absence of dust. Although the fractional abun- 
dance of H2 that forms in this way is small, it is sufficient 
to allow for effective cooling and the formation of stars 



(Haiman et al. 



L996e 



199E ; Abel et al 



2000) 



Tegmark et al 



1997 



Bromm et al. 



As soon as massive stars form, however, they imme- 
diately begin to photoionize and photodissociate this H2. 
Photoionization requires photons with energies greater than 



15.4 eV, which are strongly absorbed by neutral hydrogen, 
and is only of importance within H 11 regions. Photodisso- 
ciation, by contrast, occurs through the absorption of pho- 
tons in the Lyman- Werner band system ( stecher & Williams 
1967), with energies in the range 11.18 - 13.6 eV. These 



photons are not strongly absorbed by neutral hydro gen and 
can readily e s cape into the intergal actic medium ( jOmukai 
" |i99s| ; |Glover fc Brand| pOOll ) . 



Nishi 



Initially, many of these photons will be absorbed by in- 
tergalactic H2, but its abund ance is small and it is rapidly 



photodissociated (see section 4.4 below). Consequently, the 



onset of star formation is soon followed by the appearance 
of an ultraviolet background radiation field. This ultraviolet 
background acts to suppress further star formation by pho- 
todissociating H 2 within newly-forming protogalaxies. The 



authors (Ciardi et al. 


2000a 


200C ; 


Machacek et al 


2001 



Haiman et al 



1996b 



1997 



2001). In particular, Haiman et al.^ 
( pOOOj ) study the coupled problem of the evolution of the 
ultraviolet background and its feedback on the global star 
formation rate using a simple gala xy formation model base d 
on the Press-Schechter formalism (Press & Schechter 1974). 
They find that cooling (and hence star formation) within 
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small protogalaxies is completely suppressed prior to cosmo- 
logical reionization. Taken at face value, their results suggest 
that star formation within small, H^-cooled protogalaxies is 
a transient phenomenon, with little impact on later stages 
of galaxy formation. 

These conclusions, however, rest on the assumption that 
the only free electrons present in the protogalactic gas come 
from the small residual fraction remaining after cosmologi- 
cal recombination. This is important, as free electrons (and 
protons) catalyze H2 formation, as we can see from equa- 
tions |l| to |j. If the free electron abundance were significantly 
higher, then the H2 formation rate would also be higher, 
offsetting the effects of photodissociation. At the very least, 
this would delay the suppression of star formation, and in 
principle could entirely negate it. 

It is therefore important to determine whether there 
is any way in which an enhanced level of ionization could 
be produced. One possible source is the residual ionization 
that would remain after the recombination of H n regions 
produced by an earlier generation of stars. This has recently 
been studied by Ricotti et al. ( 2002a, who find that it can 
be an effective source of H2 and can dramatically reduce the 
effectiveness of photodissociation feedback. However, their 
conclusions are still somewhat uncertain, both becau se their 
simu lations are under-resolved (see figures 5-7 in Ricotti 
et al. 2)02a) and because they neglect the effects of super- 



novae, which would act to disperse the gas and significantly 
lengthen the recombination timescale, thereby delaying the 
formation of H2. 

An alternative possibility is ionization by a high red- 
shift X-ray background. At X-ray energies, the optical depth 
of the intergalactic medium (IGM) is small and any X- 
ray sources present will naturally generate an X-ray back- 
ground. Moreover, X-rays can penetrate to large depths 
within newly-formed protogalaxies, allowing them to raise 
the fractional ionization throughout the gas. 

The potential importance of such a background was 
first highlighted by Haiman et al. (2000), with high redshift 



quasars suggested as a possible source. Using a very sim- 
ple model for a quasar-produced X-ray background, they 
showed that if quasars contribute more than 10% of the UV 
background, then the ionization produced by the associated 
X-ray background is sufficient to negate the effects of UV 
photodissociation. Indeed, they found evidence that such a 
background could actually promote cooling within the dense 
gas in the centres of protogalaxies. This scenario has not yet 
been firmly ruled out, but observational evidence suggests 
that quasars are unlikely t o be present in suffi cient number 
at high redshift ( ^ei||l995| ; [Haiman et~al\ |l999| ). 

Quasars, however, are not the only potential source 
of X-rays. Star formation also leads to the production 
of X-rays, pri marily through the form ation of massive X- 
ray binaries (Helfand & Moran 2001), but also through 



bremsstrahlung and inverse Compton emission from super- 
nova remnants. These sources generat e only a small fraction 



of th e present-day X-ray background ( Nfatarajan fc Almaini 



200C), but may become dominant at high redshifts. 

In this paper, we examine the effects of the X-rays pro- 
duced by these sources, both on the cooling of gas within 
virialized protogalaxies and also on the chemical and ther- 
mal evolution of the IGM. The outline of the paper is as 
follows. In section 0, we discuss the sources responsible for 



producing the UV and X-ray backgrounds, and show how 
the build-up of these backgrounds can be computed. In sec- 
tion ||, we outline the method used to study the effect of this 
radiation on the primordial gas, and in section ^ apply it for 
a number of different X-ray source models. We present our 
conclusions in section n. 



2 THE UV AND X-RAY BACKGROUNDS 

For an observer at redshift zo, we can write the mean spe- 
cific intensity of the radiation bac kground at an observed 
frequency vq as (Madau et al. 199E) 



ii \ if/ \ -T{v a ,7. Q ,z) (1 + z o) 
J^,z ) = - e{u,z)e TTTzjz 



_d[ 

(1 + z) 3 di 



dz. 



(5) 



where v — Vo(l + z)/(l + zo), e(u,z) is the proper space- 
averaged volume emissivity, t(po,zq,z) is the optical depth 
at frequency due to material along the line of sight from 
redshift zo to z and Al/ dz is the cosmological line element. 
To solve this equation, we need to know how the emissivity 
and opacity evolve with redshift. 



2.1 Emissivity 

For simplicity, we write the space-averaged emissivity in 
terms of the global star formation rate (SFR) as 



s(u,z) = L(v,z)M* ergs 1 Mpc 3 Hz" 



(6) 



where L(y, z) is the luminosity density (i.e. the luminosity 
per unit frequency) per unit star formation rate (in solar 
masses per year), and M, is the global star formation rate, 
with units M Q yr" 1 Mpc -3 . 

In principle, L(u, z) may be a complicated function of 
frequency and redshift. In practice, however, Lyman- Werner 
band emission is dominated by massive, short-lived OB stars 
and declines rapidly once star formation comes to an end 
see, for example, the instantaneous starburst models of 



therer et al. 



1999 



Lei- 



in which the Lyman- Werner flux declines 
by an order of magnitude within 4-5 Myr) . Similarly, X-ray 
emission is dominated by short-lived sources such as massive 
X-ray binaries and supernova remnants which are end prod- 
ucts of the same massive stars. Both kinds of emission are 
therefore strongly correlated with the star formation rate; 
to a first approximation, we can assume that they are di- 
rectly proportional to it, and that any redshift dependence 
of L(v, z) can be neglected. With this simplification, deter- 
mining the emissivity breaks down into two independent 
problems: determining the global star formation rate as a 
function of redshift, and determining the luminosity density 
as a function of the star formation rate. 



2.1.1 The star formation rate 

Although we have observational constraints on the star for- 
mation rate up to z ~ 5, we have no direct constraints 
(and few indirect ones) at higher redshift. Consequently, 
any model of high redshift star formation must inevitably 
be highly theoretical. Moreover, this lack of constraints mo- 
tivates us to choose as simple a model as possible; more 
complicated (and realistic) models can always be considered 
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Haiman 



once our observational knowledge improves. A Roo d exam- 
ple o f this k ind of simple model is the one used by 
et al. (|j00(| ); we adopt the same model here. 

We assume that star formation proceeds primarily 
through starbursts, of duration t on years, that are triggered 
when galaxies form. During the starburst, the star formation 
rate is assumed to be constant. The global star formation 
rate in this model is given by 



e /bAp ga l 



M 



y r 



1 Mpc" 



(7) 



where e is the star formation efficiency, /b is the baryon 
fraction (ie the ratio of baryons to dark matter), and where 
Apgai is the cosmological density of matter in newly-formed 
galaxies (with units of M Q Mpc -3 ). We assume that the 
value of /b in the protogalaxies is the same as in the IGM, 
or in other words that /b = Qf,/fi m . For the cosmological 
model adopted in section ^, this corresponds to /b = 0.124. 

We further assume that the rate of change of Ap ga i is 
approximately the same as the rate of change of F[z, T cr i t ), 
the total fraction of matter in halos with virial temperatures 
greater than a critical temperature T cr it, or in other words 
that 



dz 



(8) 



where p m = 2.8 x 10 n n m (l + z) 3 M© Mpc -3 is the cosmo- 
logical matter density. Here, T cr it represents the minimum 
virial temperature required for efficient cooling; to a first ap- 
proximation, halos with T v ir < Tcrit are unable to cool, while 
those with T v i r > T cr i t cool rapidly and can form stars. Var- 
ious different de finitions of T cr i t a re in use in the literature 
(see, for example [Rees fc Ostrikcr 



Haiman et al 



1977 



Tegmark et al 



1997 



2000[ ); we discuss our particular choice in a 



later section 

Making this approximation is equivalent to assuming 
that the growth in F(z, T cr i t ) is dominated by the formation 
of new halos with T v i r ^ T cr it (either by monolithic collapse 
or by the merger of smaller objects) rather than by the ac- 
cretion of matter by existing halos with T v j r > T cr i t . This 
is justified at high redshift when such halos are rare and 
F(z, T cr it) is dominated by objects near Tcrit, but becomes 
less accurate at lower redshifts. For this reason (and others, 
to be discussed later) we do not attempt to simulate the 
evolution of the background below z = 10. 

To solve equation H we need to know the value of T C rit ■ In 
gene ral, th is will depend both on redshift (see, eg Tegmark 
et al. 1! 97) and on the intensities of the UV and X-ray back- 
grounds. Indeed, understanding the evolution of T cr it with 
redshift is one of the main goals of this paper. Our procedure 
for determining T cr it is discussed at length in section [| for 
now, we assume that it is known. In this case, we can calcu- 



late F (z, Tcri t) using the Press- Schechter formalism (Lacey 
& Cole j|l993[ ): 



F(z, Merit) = erfc 



S c (z) 



%/2<r(Mcrit) 



(9) 



where 5 c (z) is the critical density threshold for collapse, 
a(M) is the square root of the variance of the cosmologi- 
cal density field, as smoothed on a mass scale M, and where 
Merit is the mass of a protogalaxy with virial temperature 
Tcrit- Although both S c and <r(M) depend upon the choice of 



cosmological model, their behaviour is well-known, and the 
problem of determining F(z, M cr it) reduces to the relatively 
simple one of relating Af C rit to T cr it- 

To do this, we need to know the protogalactic density 
profile. Hydrodynamical simulations sug gest that it is ap- 
proximately isothermal (Abel et al. 2000| ), but representing 
it as a singular isothermal sphere is physically unrealistic 
due t o the latter's in finite central density. Accordingly, we 
follow 



Haiman et al 



and represen t it as a truncated isother- 
mal sphere ( |llicv & Shapiro] 2001). With this choice, we find 
that 



Merit = 1.5 xlO 7 (l+^r^n-^/^Me-ao) 

To determine Ap ga i(z), we must integrate equation ^ 
over a redshift interval Az(t a n), corresponding to the dura- 
tion of the starburst; hence, 



Ap g ai(2) 



: + Az(t on ) 



dF 

pm{z)——(z, Tcrit) dz. 

dz 



(11) 



This simple model has a number of shortcomings. For 
instance, it assumes that the star formation efficiency e and 
starburst duration ion are both constant, independent of red- 
shift or galaxy mass. It also assumes that each galaxy forms 
its stars in a single starburst and thus ignores the effects 
of continuous star formation and of subsequent, merger- 
triggered starbursts. [] Nevertheless, it has the virtue of sim- 
plicity, and is a good point from which to start our exami- 
nation of the effects of the X-ray background. 



2.1.2 The UV luminosity density 

The ultraviolet flux of a star-forming galaxy is dominated 
by emission from young, massive O and B-type stars. These 
are short-lived, with the most massive having lifetimes of 
only a few Myr, and thus the ultraviolet luminosity density 
is closely correlated with the star formation rate. Its value 
depends upon the spectral properties of the newly-formed 
stellar population, and thus on their initial mass function 
(IMF), metallicity and age. 



In a recent paper, Schaerei (2002) presents values for 



the photon flux in the Lyman- Werner bands calculated for 
a number of different metal-free stellar populations. If we 
assume that the spectrum within the bands is flat (a reason- 
able approximation), then we can convert this photon flux 
into a luminosity density. For a Salpeter IMF with minimum 
mass M m i n — 1 Mr and m aximum mass M m ax = 100 Mq 
(model A in |5chaererj |5o§) , we find that 



L v = 1.1 x 10 



28 ergs" 1 Hz -1 



M 



y r 



(12) 



Reducing M m i n to the more conventional value of 0.1 Mq 
reduces this luminosity to 



L v = 4.3 x 10 2 ' ere 



^Hz" 1 (Moyr- 1 ) \ 



(13) 



as we form a greater number of low mass stars that do not 
contribute significantly to the Lyman- Werner flux. 

Both of these results assume that the dissociative flux 



1 Although both of these can be mimicked to some extent by 
choosing a large value for t on . 
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has stabilized at its equilibrium value and is therefore pro- 
portional to the star formation rate. This equilibrium is typ- 
ically established after only 2-3 Myr, so this is generally a 
good approximation, even for starbursts of relatively short 
duration. 

The above figures are appropriate so long as we are deal- 
ing with stars formed out of entirely metal-free gas. Such 
stars are somewhat unusual, however, as the absence of car- 
bon means that they are unable to generate energy via the 
CNO cycle, which otherwise would dominate energy pro- 
duction in stars of mass M > 1.1 Mq. As a result, metal-free 



stars ar e hotter th a n their metal-enriched counterparts ( Ezci 



Camsron 



Shuf 



2000 



11971 ;_ Sassisi & Castellan: 



1993 



Tumlinson 



Sojazzi et al. 200C) and have harder spectra. A 



surprising consequence of this fact is that a low-metallicity 
stellar population will produce a larger dissociative flux than 
a metal-free population - the lower effective temperatures 
move the peak in the thermal emission from the most mas- 
sive stars closer to the Lyman- Werner band, causing the 
ionizing flux to fall but the dissociativ e flux to r ise. 



We can use the data presented in Bchaerer (2002) for a 
stellar population with Z — 0.02Zq to examine the differ- 
ence that this effect makes to the Lyman- Werner flux. Using 
the same IMF as in equation [H| we find that 



L v = 1.8 x 10 28 erg s _1 Hz -1 (M Q yr _1 ) " 



(14) 



Thus, raising the metallicity increases the Lyman- We rner 
flux, but only by about 60-70%; as we will see in section 4.1 
this has little effect on the evolution of T cr it. 

Clearly, there are many possible models other than 
those considered here. Indeed, there is growing evidence that 
the IMF of p opulation III stars is strongly biased towards 



high masse s ( Larson 199ij ; Abel et al. 



199£ 



2000 



Bromm et al. 



200S). However, this remains uncertain, and in this 



paper we have chosen to err on the side of caution and as- 
sume that the high-redshift IMF is similar to that at the 
present day. 



2.1.3 The X-ray luminosity density 



In a recent study, Helfand & Morar (2001) collate data on 



a number of local starburst galaxies and compare the 2 - 
10 keV X-ray fluxes measured by ASCA with the 8 - 1000 
/j,m infrared fluxes measured by IRAS. They find that a clear 
correlation exists, with the total X-ray and infrared fluxes 
related by 



Fx ~ 10" 4 Fi R . 



(15) 



Simil a r corr elati ons have previously be en reported by David 



et al. 



( |1992{ ) and [Rephaeli etal\ ( |l995[ ) for X-rays in the 0.5 
- 4.5 keV and 2-30 keV energy bands respectively. 

Theoretically, we would expect such a correlation, with 
both the X-ray and infrared emission tracing the underlying 
star formation rate. For the X-rays, this occurs because the 
emission is dominated by massive X-ray binaries (MXRBs): 
binary systems consisting of a massive OB star accreting 
onto a compact companion (a neutron star or black hole). 
X-ray emission from such systems generally switches on a 
few million years after the formation of the compact object, 
and the lifetime of the emitting phase is short, typically 



of the order of (2 - 5) x 10 4 yr ( Meurs fc van den Heuvefl 



198£). These short timescales ti e the emission closely to the 
underlying star formation rate (Ghosh & White [2001 ). The 
far-infrared flux, on the other hand, tracks star formation 
far more directly, being dominated by emission from dust 
heated by short-lived, massive stars. 

To use this observed correlation to determine the X-ray 
luminosity of a star-forming galaxy as a function of its star 
fo rmation rate, we use the resu lt from the starburst models 
of Lcitherer & Heckman (199E) that 



1.5 x 10 1 



SFR 



L©, 



(16) 



1 M yr- 

t ogether with equation lj| to write the X-ray luminosity as 
flHelfand & Moran||20uT) 



Lx = 6 x Iff 3 



SFR 



1M 



© yr" 



erg s 



(17) 



We then assume that this relationship between X-ray lumi- 
nosity and star formation rate remains valid as we move to 
higher redshifts. Note that the same need not be true for the 
relationship between X-ray luminosity and infrared luminos- 
ity, or infrared luminosity and star formation rate: although 
they serve to establish the correlation between X-ray lumi- 
nosity and star formation rate at z = 0, the assumption that 
this correlation remains valid at higher redshift does not im- 
ply that these other correlations also remain valid. Indeed, 
we would expect the infrared luminosity of dust-free proto- 
galaxies to be very much lower than would be predicted by 
equation 

Evidence that the correlation between X-ray luminosity 
and star formation rate does indeed remain valid at high red- 
shift is provided by the recent stacking analysis of individu- 
ally undetect ed Lyman break g alaxies in the Chandra Deep 
Field-North ( [Brandt et a7|[>00l| ). This analysis finds that the 
average rest frame luminosity of the Lyman break galaxies 
in the 2-8 keV energy band is Lx = 3.2 x 10 41 erg s _1 . 
Assuming a ty pical star formation rate of 50 Mq yr -1 for 



these galaxies (Shapley et al. 2001), this corresponds to an 
X-ray luminosity of 



L x = 6.4 x 10 d 



SFR 



1M 



© y r " 



erg s 



(18) 



consistent with the value derived above. f\ Although far from 
conclusive, this result suggests that we can extrapolate the 
locally observed correlation to at least as far as z ~ 4. 

Compa ring our de termination of Lx with a recent cal- 
culation by Oh ( 2001 ), we find a difference of a factor of 
ten in our results. Some of this disagreement is due to the 
diff e rence in X-ray energy bands considered (0.2 - lOkeV in 
Oh (2001), compared to 2 - 10 keV here), but some must 
surely be due to intrinsic scatter in the observational data, 
suggesting that equation [lj] should properly be regarded as 
an order of magnitude estimate of the true X-ray luminosity. 



2 Note that X-ray luminosity of equation |l8] is measured in a 
slightly different energy band from that of equation so the 
agreement between the two values is not quite as good as may at 
first appear. Nevertheless, the necessary correction is small, and 
the values agree to within 50%. 
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Given equation 0for the X-ray luminosity as a function 
of the star formation rate, we calculate the X-ray luminosity 
density by assuming a template spectrum of power-law form 



L„ — L 
where hvo 

"V 2 



L x = 



•(?)'■ 

= 1 keV, and requiring that 
L„ &v, 



(19) 



(20) 



where hv\ — 2 keV a nd hv 2 = 10 keV. 

find that a weighted average of 



Rephaeli et al. (199 



the galaxies in their sample gives a value for the spectral 
index of a = 1.5 ± 0.3. Adopting this value and solving for 
Lq, we find that 



L = 3.4 x 10 22 ergs -1 Hz -1 (M 



(21) 



Altering a changes Lq, but never by m ore than 50% for val 



Rephaeli et al. measurement. Clearly, 



ues consistent with the 
individual galaxies may have spectra that differ markedly 
from this simple template, but it should be a reasonable ap- 
proximation when averaging over a large number of galaxies. 



2.1-4 Alternative X-ray models 

The above model is simple, and empirically motivated, but 
does assume that the X-ray emission of high-redshift star- 
forming galaxies is very similar to that of starbursts observed 
locally. This is a reasonable assumption in the absence of 
evidence to the contrary, and is probably valid as long as 
massive X-ray binaries continue to dominate the galactic 
X-ray emission. However, it is quite possible that at high 
redshift some other type of source will come to dominate 
the emission, particularly if the number of bin ary systems 



is sm all, as is suggested by recent simulations (Abel et al. 



200C). It is therefore prudent to consider the effects of other 
potential sources of X-rays. 

The obvious candidates are supernova remnants (SNR); 
next to X-ra y binaries, they are the m ost significant galac- 
tic sources (Helfand & Moran 2001). They can emit X- 
rays through a variety of different emission mechanisms, 
but at high redshift the most significant will be thermal 
bremsstrahlung emission and non-thermal inverse Compton 
emission. 

Thermal bremsstrahlung is produced by the hot gas 
within the SNR. Detailed modeling proper ly requires a hy - 
drodynamical treatment (see, for example, Chevalier 1999), 
but for our purposes a simple parameterization suffices. If 
we assume that the hot gas has a single characteristic tem- 
perature T x , then we can write the luminosity density per 
unit star formation rate as 

L v = L b exp (=g) ergs^Hz- 1 (M yr" 1 )" 1 .(22) 

where hvo = 1 keV and where Lb is constant. Moreover, we 
can write the total X-ray luminosity as 



Ltot = 3.2 x 10 43 / x £ B iiV sn ergs" 1 (M Q yr 



(23) 



where £51 is the typical supernova explosion energy (in units 
of 10 51 ergs), / x is the fraction of this energy radiated as 
bremsstrahlung, and where N sn is the number of supernovae 
that explode per solar mass of stars formed. The value of N Bn 



depends on the IMF; for the standard Salpeter IMF adopted 
previously, N Bn — 0.0075 Mg 1 . Finally, since 



Ltot = I Lv Av, 
10 



(24) 



we can fix the value of Lb; we find that 
L b = 1.3 x 10 26 / x L 5 i A^sn erg s" 1 Hz" 1 (M Q yr" 1 ) _1 . (25) 
For typical supernova p arameters (L51 = 1 and an am- 



Helfand & Moran 



(2001) find 



bient density n = 1 cm ), 
that a fraction f K — 2 x 10 -4 of the explosion energy is ra- 
diated, at a characteristic temperature fcT x = IkeV. On the 
other hand, the higher mean density at high redshift, to- 
gether with the c omparative weak ness of outflows from low 
metallicity stars (Kudritzki 2000) suggest that the typical 
ambient density may be very much higher. In particular, if 
it is as high as n ~ 10 7 cm~ 3 , then a supernova remnant 
will radiate its energy e xtremely rapidly, bef ore the ejecta 
have time to thermalize (Terlevich et al. 1992). In this case, 
the fraction of energy radiated as bremsstrahlung is very 
much higher (/ x ~ 0.01), as is the characteristic tempera- 
ture (fcT x ~ 30 keV). We examine both of these models in 
section ^| with the understanding that the true picture lies 
somewhere in between. 

In addition to this thermal emission, supernova rem- 
nants also produce non-thermal X-rays. These are generated 
as the relativistic electrons produced by the SNR gradu- 
ally lose energy through synchrotron radiation, non-thermal 
bremsstrahlung and the inverse Compton scattering of pho- 
tons from the cosmic microwave back ground. A t high red- 
shift, the latter is likely to dominate (Ob 2001). The spec- 
trum of the resulting emission depends upon the energy 
spectrum of the relativistic electrons, but at the energies 
of interest is well represented by a power law: L„ oc v^ 1 . 
The intensity of the emission depends upon the fraction of 
the supernova energy transferred to the electrons; this is not 
well constrained, with estimates ranging from 0.1% to 10%. 
Accordingly, we model the emission as 

L u = 7.7 x 10 23 /. erg s' 1 Hz" 1 (Mq yr" 1 ) ~ X , (26) 

where hvo — 1 keV, f e is the fraction of energy deposited in 
the electrons, and where we have assumed that L51 = 1 and 
A sn = 0.0075 Mq 1 as in the thermal bremsstrahlung case. 
This expression assumes a high-energy cutoff for the X-ray 
spectrum at 10 keV but is only logarithmically dependent 
on the value of this cutoff. In section ^, we examine results 
for models with f c = 10~ 3 and f c — 0.1, which bracket the 
range of plausible values. 

Much more detail on high redshift inverse Compton 
emission, includi ng a discu ssion of potential observational 
tests, is given in Oh (2001). 



2.2 Opacity 

The opacity t(i/q, zq, z) can be separated into two distinct 
components - absorption by dust and gas within the emit- 
ting protogalaxy, which we term intrinsic absorption, and 
absorption by gas along the line of sight through the IGM. 

Intrinsic absorption is difficult to model with any de- 
gree of accuracy as it depends upon a number of variables - 
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the size and shape of the galaxy, its ionization state, the po- 
sition of the sources within it, the dust content etc. Rather 
than attempt to model these in detail - a significant under- 
taking in itself - we instead adopt a highly approximate rep- 
resentation. We assume that the emitted X-rays are attenu- 
ated by absorption by a neutral hydrogen column density of 
Nb. = 10 21 cm -2 plus an associated neutral helium column 
density 7Vh c = 0.08A?k- This absorption is assumed to be 



1 + Z 



the same for all sources. These values are chosen because 
they are representative of the column densities of the proto- 
galaxies studied in this paper (which presumably contain the 
bulk of the X-ray sources). Reducing JVh (as would be ap- 
propriate if much of the surrounding Hi were photoionized 
and/or dispersed by the progenitors of the X-ray sources) 
has little or no effect on T cr i t , as gas in the IGM and within 
the protogalaxy itself quickly come to dominate the total 
absorption. Increasing JVh, on the other hand, has more sig- 
nificant effects: an order of magnitude increase in TVh pro- 
duces similar results to an order of magnitude decrease in 
the strength of the X-ray background, which, as we will see 
in section [|, is generally sufficient to render X-ray feedback 
ineffective. Consequently, we will overestimate the effect of 
the background if the bulk of the X-ray sources reside within 
massive galaxies. 

Turning to the Lyman- Werner bands, we note that in- 
trinsic absorption will generally be negligible within small 
proto galaxies, as their H2 content is rapidly photodissoci- 
ated ( Glover fc Brancj 2001). Moreover, we also assume that 
the effects of dust absorption are negligible. In galaxies of 
primordial composition, this is obviously true; in metal-poor 
galaxies, it should also be a good approximation, as very 
large column densities are required for significant dust ob- 
scuration (for instance, TVh ^ 10 23 cm -2 for Z — W~ 2 Zq 
gas, if the dust-to-gas ratio is similar to that in the Milky 
Way). Again, these assumptions break down if the majority 
of sources are to be found in massive, metal-rich galaxies, 
but we expect such galaxies to be extremely rare at the red- 
shifts of interest in this paper. 

Compared to intrinsic absorption, the effects of absorp- 
tion due to the IGM are much simpler to treat, particularly if 
we can assume that the bulk of the gas remains at an approx- 
imately uniform density. This assumption proves reasonable 
at high redshift for photons with mean free paths much 
greater than the typical clumping scale, as is the case for 
both Lyman- Werner band photons and X-rays. Our treat- 
ment of IGM absorption is discussed in the following sec- 
tions. 



2.2.1 Lyman- Werner band absorption 

The continuum op acity of metal-free gas is very small 
(Lenzuni et al. 1991) and for our purposes can be neglected. 
Consequently, the only significant sources of opacity encoun- 
tered by Lyman- Werner photons are absorption by the Ly- 
man series lines of neutral hydrogen, and by the Lyman- 
Werner lines of molecular hydrogen. The Lyman series lines 
have the effect of absorbing any Lyman- Werner photons of 
the same frequency, and reprocessing them to Lyman-a pho- 
tons plus associated softer photons. As Lyman-a lies outside 
of the Lyman- Werner band, the net effect is to block from 
view any sources at redshifts higher than some maximum, 
given by 



1 



max 



■ Zq 



(27) 



where is the frequency of the appropriate Lyman series 
line and fo and zo are the observed frequency and redshift. 
Clearly, the size of z max depends upon the distance between 
vo and vh , and thus more sources are seen at frequencies that 
are a long way from a line. As a result, the spectrum develops 



a character istic 'sawtooth' shape (see figure 1 in Haiman 
2000), with the effect becoming more pronounced as 



et al 



one nears the Lyman limit. 

Absorption by molecular hydrogen is rather more com- 
plicated, due to the large number of Lyman- Werner lines 
that contribute to the opacity. If we approximate the lines 
as delta functions then an individual line produces an opac- 
ity 



" /osc ,i -V 



H{ Zi ) 



(28) 



(29) 



where / OS c,i and Xi are the oscillator strength and wave- 
length of the transition, nH 2 ,i( z i) is the number density of 
H2 molecules in the level giving rise to the line at Z4, the 
redshift of absorption, and H(zi) is the Hubble constant at 
Zi. The value of Zi is given by 

1 + Zi _ Vi 

1 + zo uo ' 

where Vi = c/Xi. 

If we assume that all of the photons that are absorbed in 
the lines are permanently removed from the Lyman- Werner 
band, then the total opacity t(vq, zq, z) is simply given by 
the sum over all lines with z < Zi < z max : 

ne ' 2 r \ n K 2 , t (zi) , . 

T{V ,Z0,z)=y fosc.iXi . (30) 

^ — ' m,„c. H { Zi\ 



This sum potentially involves a very large number of lines, 
but can be greatly simplified by assuming that all of the 
H2 is to be found in its ortho or para ground state; at the 
redshifts of interest, the population of excited states will be 
negligible. 

In deriving this expression, we have assumed that ev- 
ery absorption permanently removes a Lyman- Werner band 
photon. This is not entirely correct. On average, only 15% 
of absorp tions are followed by ph otodissociation of the H2 



molecule ( Draine fc Bertoldi 1996 ); the rest of the time, the 



molecule decays back to a bound sta te, emitting a photon 



In their treatment of this problem, Haiman et al. (200C) 
assumed that the excited H2 molecule would always decay 
directly back into the original state, and would thus emit 
a photon of the same energy as the one initially absorbed. 
In fact, this is not correct (T. Abel, private communica- 
tion); most decays occur initially to highly excited vibra- 
tional states, producing photons redwards of the Lyman- 
Werner bands. Only a small fraction of decays (about 5%) 
take place directly into the original state, while a slightly 
larger fraction (about 15%) produce phot ons that lie else- 



wher e in the Lyman- Werner band system (Glover & Brand 
2001). We do not include the effect of these photons, how- 
ever; an accurate treatment would be quite complicated and 
is almost certainly unnecessary - as we shall see in section ^j, 
H2 in the IGM is rapidly destroyed by the growing Lyman- 
Werner background and is completely negligible by the time 
that negative feedback begins. 
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2.2.2 X-ray absorption 

At X-ray energies, the opacity of the intergalactic gas is 
dominated by the ionization of neutral hydrogen and helium; 
prior to recombination, the He + abundance is small and can 
be neglected. The X-ray opacity can thus be written as 



Tcrit ( z); th e de tails of this simulation are outlined in sec- 



"(fo, «o, 



[(Til(v) n H + (Tn c (u) 7lHe 



dz 



dz. 



where 



v = Vq 



(1 



(l + *o)' 



(31) 



(32) 



and where an(u) and o~Ke{y) are the absorption cross- 
sections of neutral hydrogen and helium respectively, with 
nn(z) and ?iHe(z) being the corresponding number densi- 
ties. As long as the fractional ionization of the IGM remains 
small (i.e. a few percent or less), the ratio between wh c to 
n.H can be accurately approximated by its primordial value 

— - Y (33) 

where Y = 0.247 is the helium mass fraction, and we can 
write equation |3l] purely in terms of nu as 

Y 



t(v , z , 



£>~H + 



(4-4^) 



0~Hc 



nu—dz. 

dz 



(34) 



This integral is readily computable by means of numerical 
integration. 



3 METHOD 

In the previous section, we showed that, given a simple star 
formation model, it is relatively easy to calculate the evolu- 
tion of the Lyman- Werner and X-ray backgrounds. Two of 
the parameters of our star formation model - the star for- 
mation efficiency e and starburst lifetime t on - we treat as 
free parameters ( although they can be constrained to some 



extent - see, e.g. Haiman & Loet 1997). The remaining pa- 
rameter, T cr it, the temperature at which efficient cooling be- 
comes possible, is determined by the strength of the back- 
grounds themselves. This clearly presents us with a problem: 
the evolution of T cr i t is coupled to that of the backgrounds, 
and to know one we must first know the other. 

Fortunately, this difficulty is easily avoided. We know 
that at high redshift the number of protogalaxies, and hence 
the star formation rate, must be very small. Consequently, 
there must be some redshift above which the external radi- 
ation field will become too weak to affect galaxy formation. 
The precise redshift at which this occurs is model dependent, 
but for the models examined in this paper we typically find 
that radiative feedback is negligible above 2 = 40. By choos- 
ing an initial redshift Zi = 50, therefore, we can be sure that 
in our initial simulation the background radiation will have 
no effect. 

Given this starting point, we next proceed incrementally 
to lower redshifts via the following procedure: 

(i) Given T C iit(Zi), we calculate the background radiation 
field at z — Zi — Az, assuming that T CI it(z) = T CI i t (zi). 

(ii) Using the computed background, we simulate the 
chemical and thermal evolution of a protogalaxy with T v i r = 



tions 3.1 



to 3.7 



below. The main aim of this simulation is to 
determine whether the protogalactic gas can cool efficiently. 

(iii) If the protogalactic gas cools, then our assumed value 
of Tcrit is correct; we store this result, and return to step one 
to proceed to the next redshift. If the gas fails to cool, we 
continue to step four. 

(iv) We increment our assumed value of T C rit (z) by a small 
amount AT, and recalculate the background radiation field. 
We assume that T cr i t varies linearly over Az. Given the new 
background, we return to step two. 

Provided that Az and AT are both small, the error in 
Tcrit (z) will also be small; this is particularly the case once 
emission from larger protogalaxies (which cool via Lyman-Q 
radiation) begins to dominate the background. 

This approach reduces the coupled problem to the sim- 
pler one of determining whether a protogalaxy with virial 
temperature T v i r and formation redshift z; will cool when 
exposed to a particular background radiation field. To an- 
swer this question, we need to be able to model the thermal 
and chemical evolution of the protogalaxy. Our approach to 
this problem is outlined in the following sections. 



3.1 Computing the evolution of protogalactic gas 

Ideally, we would like to use a high-resolution hydrodynam- 
ical simulation to follow the thermal and chemical evolution 
of the protogalactic gas (see, e.g. 



Machacek et al. 



2001). Un- 



fortunately, including the effects of radiative transfer, par- 
ticularly of photons in the Lyman- Werner bands, into such 
a simulation is not currently feasible. We are thus forced 
to approximate. In choosing an appropriate approximation, 
we are also motivated by the desire to minimize the com- 
putational requirements of our simulations, so that we can 
explore the effects of a variety of different source models. 
We make three main approximations: 

(i) We assume spherical symmetry. This is a reasonable 
approximation for the first generation of protog alaxies, b ut 



clear ly is incorrect in detail (see, e.g. figure 2 of Abel et al 
200C|). 



(ii) We assume that the protogalactic gas is static, at least 
on the timescale of the simulation. This allows us to ignore 
the hydrodynamical evolution of the gas, and also substan- 
tially simplifies the treatment of radiative transfer. This as- 
sumption clearly breaks down once the gas begins to cool 
strongly and loses its pressure support, but as we are only 
interested in the evolution up to this point, this is not a 
significant problem. 

(iii) We assume that all of the H2 molecules remain in 
the rotational and vibrational ground state, in either ortho 
or para form. This simplification allows us to evolve the 
chemistry and radiative transfer on the timescale on which 
the total H2 abundance changes (typically 10 11 - 10 12 s) 
rather than that on which the individual level populations 
change (10 6 - 10 10 s). It also simplifies our treatment of the 
radiative transfe r. We discuss this approximation in more 
detail in section 3.4.2. 



Together, these approximations allow us to solve for the 
chemical and thermal evolution of a model protogalaxy in a 
matter of minutes on a fast desktop computer. This allows us 
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to study the redshift evolution of T cr i t at high resolution in 
both temperature and redshift and for a number of different 
X-ray source models. 

However, this approach has an obvious drawback - we 
cannot be sure that our approximations give a fair repre- 
sentation of the real protogalaxy. Of particular concern is 
the neglect of the hydrodynamic evolution of the gas, and 
the consequent error in the density profile. This is poten- 
tially significant because the H2 cooling rate, along with 
many of the chemical reaction rates, scales as the square 
of the density. Small errors in the density can thus lead 
to larger ones in the computed temperature. These con- 
cerns are mitigated to some extent, however, by the close 
agreement between the results of detailed numerical simu- 
lations and previous semi-analytic treatments. For instance, 
the values of T cr it obtained from the smoothed - parti cle hy- 
drodynamics simulations ofJ^Jki^^^oujdunan (200C) agree 
well with the results of Tegmark et al. ( 1997 ), despite the 
highly approximate uniform density profile adopted by th e 
latter group. Similarly, the results of Machacek et al. (2001), 
obtained with a three-dimensional adaptive mesh hydrody 



namic al code broadly agree with those of Haiman et al. 
(2000), who use a static model similar to that presented 
here. Together, these results suggest that T cr it is insensitive 
to the precise details of the density profile, but clearly this 
remains an area of concern. 

Our computational method can be broken down into 
three main stages - initialization of the density profile and 
the chemical abundances, computation of the thermal and 
chemical evolution of the gas, and termination of the simu- 
lation at a suitable point. These are described below. 



3.2 The model protogalaxy 

The protogalactic density profile is modeled as a truncated 
isothermal sphere, with central overdensity S = 1.796 x 10 4 
and truncation radius 



r t = 4.62 x W 



The virial temperature and redshift of formation of the pro- 
togalaxy completely specify its density profile. 

We subdivide this profile into N g spherical shells of uni- 
form thickness and compute the mean density within each 
shell. We have run a number of test simulations with differ- 
ent values of N g , and find that setting N g — 100 provides 
sufficient spatial resolution to accurately determine T cr i t . 

We assume that the initial chemical composition of the 
protogalactic gas is the same as that of the intergalactic 
medium. At our ini tial redshift (z — 50), we take this from 
Standi et al. ( 199S ) . At lower redshifts, the chemical evo- 



lution of the IGM is influenced by the Lyman- Werner and 
X-ray backgrounds. We therefore calculate the intergalactic 
abundances explicitly, using the chemical model outlined in 
section 3.3, by solving the chemical rate equations: 



drii 



Ci{n h T)-Di(n h T)ni 



(36) 



where d and D; are source and sink terms for rii. At the 
same time w e also solve for the temperature of the inter- 
galactic gas (|Galli fe Paiia| ft998|) 



dT 
~~dt 



-2TH(z) + 



3kn to 



(r - A) , 



(37) 



where H(z) is the Hubble constant, n to t is the total particle 
number density and wher e F and A are the net heating and 
cooling rates (see section 3.5). 

We solve this set of ordinary d i fferen tial equations with 
the STIFBS integrator of Press et al. (1986). As the size of the 



required timestep is generally much smaller than the redshift 
interval Az that separates our individual determinations of 
Tcrit, we compute intermediate values by interpolation and 
from these determine the strength of the radiation back- 
ground and hence the photochemical rates. Although our 
main aim in following this chemistry is to determine the 
correct initial abundances for our simulations of protogalac- 
tic evolution, the results are of interest in the ir own right 
and are presented and discussed in section 4.4. 



Our treatment of the chemistry of the IGM does not in- 
clude the effects of the ionizing photons from stars (and/or 
quasars) that are ultimately responsible for the reionization 
of the intergalactic gas. This is justified at early epochs, 
as these photons are confined within small H 11 regions sur- 
rounding the luminous sources, but this simplification re- 
stricts the validity of our results to the period prior to cos- 
mological reionization. The post-reionization epoch, and the 
effect of reionization on gala xy formation, have rec eived ex- 



tensive study elsewhere (see 
erences therein). 



Loeb & Barkana 2001 



and ref- 



3.3 The chemical model 

To simulate the chemical evolution of the protogalactic gas, 
we adopt a chemical model consisting of thirty reactions 
between nine species: H, H~, H + , He, He + , He ++ , H2, H^" 
and free electrons. The reactions included in the model are 
summarized in table together with the source(s) of the 



data u sed. T his model is based in large part on that of Abe] 



et al (1997), but has been modified to improve its accuracy 
when applied to optically thick gas. Aside from a number of 
updates to the reaction coefficients in the light of new data, 
the main differences are as follows: 

(i) We include the contribution to the hydrogen ioniza- 
tion rate arising from the ionizing photons produced by He + 



recombination, in line with the discussion in chapter 2 of Os 



terbrock (1989). Although commonly a small correction to 
the total rate, this can become significant when X-ray pho- 
toionization dominates. 

(ii) To enable us to accurately determine the He + abun- 
dance, we find that we need to include the effects of charge 
transfer between He + and H (reaction 20), as this can be 
comparable to the recombination rate when the fractional 
ionization is small. For completeness we also include the 
inverse reaction (no. 21), although this is unimportant at 
T < 10 4 K. 

(iii) We include the contribution to the ionization rates of 
hydrogen and helium arising from secondary ionization by 
energetic photoe l ectro ns, based on the recent calculations of 
Dalgarno et al. (1999). The contribution of secondary ion- 



ization to the other photoionization rates is small and can 
be neglected. 

(iv) We do not include the photodissociati on of H2 b y 
photons above the Lyman limit (reaction 28 in Abel et al.), 
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No. Reaction 



Reference 



1. 


H + 


e — ► 


H+ +2e 


2. 


H+ 


+ e - 


-> H + 7 


3. 


He H 


- e — 


Hc+ + 2e 


4. 


He+ 


+ e 


-> He++ + 2e 


5. 


He+ 


+ e 


-» He + 7 


6. 


He+ 


+ + e — > He+ + 7 


7. 


H + 


e — > 


H- +7 


8. 


H- 


+ H 


^H 2 + e 


9, 


H + 


H+ - 


-H++7 


10. 


H 2 + 


+ H- 


- H 2 + H+ 


11. 


H 2 - 


f-H+ 


— H+ + H 


12. 


H 2 - 


He - 


> 2H + e 


13. 


H 2 - 


hH- 


■+ 3H 


14. 


H- 


+ e - 


^H + 2c 


15. 


H - 


+ H 


-» 2H + c 


16. 


H" 


+ H+ 


-> 2H 


17. 


H" 


+ H-+ 


^H++e 


18. 




-fc - 


* 2H 


19. 


H" 


+ He 


-> H + He + c 


20. 


He+ 


+ H 


-» He + H+ + 7 


21. 


Hc- 


-H+ 


-> He+ + H 


22. 


H + 


7 -» 


H+ +e 


23. 


He H 


-7- 


> He+ + e 


24. 


He+ 


+ 7 


-» He++ + e 


25. 


H" 


+ 7- 


-> H + e 


2(5. 


H 2 - 


1-7- 


+ H ++e 




Moseley et gq (p97C 
Poulaert et at 



27. 
28. 
29. 
.30. 



H 2 + 7 " 



11] 
11.1 



+ 7 " 
+ 7- 
H 2 + 7 - 



H + H+ - 
> H + H+ 
» 2H+ +e 



Hit 



2H 




einhardq (|197J) 
Wilms et ,a | f|200o| 



Yan et al 



Samson k Haddac ( L994 ) 



Dunn 



(196S) 



Bates fc Opik (1968) 



See text 



Table 1. A list of the reactions included in our chemical model 
of protogalactic gas. Values for the rate coeff icients f or rad iative 
cross-sections where appropriate) are given in Glover (2001). Ref- 



erences are to the primary source(s) of the data whenever possible; 
in many ca s es, we have also used analy tica l fits to this d a ta fro m 
Abel et al\ ( 1997| ), |Galli fc Palla| ( |l998[ ) or |stancil et at (1998). 



as in optically thick gas this will be negligible compared to 
the effects of H 2 photoionization (reaction 26). On the other 
hand, we do include the effects of dissociative photoioniza- 
tion (reaction 27), which becomes significant for photon en- 
ergies greater than 30 eV. 

More information about all of these points , and the 
chemical model generally, can be found in Glover (2001). 



3.4 Radiative transfer 

In order to determine rate coefficients for the photochemical 
reactions (nos. 22-30), we need to calculate the evolution of 
the mean specific intensity J„ within the protogalactic gas. 
For optically thin gas, this is simple, but the only photo- 
chemical reactions for which we can assume optically thin 
conditions are H~ photodetachment (reaction 25) and Hj 
photodissociation (reaction 28); for the rest, we must sup- 
plement our chemical model with a treatment of radiative 
transfer. 

The high degree of symmetry of the problem makes 



direct solution of the radiative transfer equation feasible, 
but the large spectral resolution required to properly treat 
the Lyman- Werner lines (Av/u ~ 10 ) makes such an ap- 
proach very time-consuming, particularly since it must be 
repeated for a very large number of timesteps. This moti- 
vates our choice of a much simpler approximate treatment 
that nevertheless should produce results of an acceptable 
accuracy. 

To begin with, we assume that the re-emission and/or 
scattering of photons within the protogalactic gas is negligi- 
ble. Although highly accurate for the X-ray photons, this 
approximation is less so for the Lyman- Werner photons, 
since many Lyman- Werner absorptions are follo wed by re- 
emission. However, as noted previously in section [2.2.1 , the 
majority of these photons are emitted at frequencies that 
do not coincide with the ground-state Lyman- Werner lines, 
and rapidly escape from the protogalaxy without being re- 
absorbed. By ignoring the small fraction of re-emitted pho- 
tons that do coincide with ground-state lines, we underesti- 
mate the H 2 photodissociation rate slightly, but this effect is 
small compared to the uncertainty in the incident spectrum. 

In addition to the ionizations produced by the incident 
X-rays (and secondary photoelectrons), we also must ac- 
count for the ionizations caused by the emission of ionizing 
photons during H and He recombination within the gas. This 
is done by use of the on-the-spot approximation, as outlined 



Osterbrock (198E). This approximation is valid as long as 
the mean free path of the ionizing photons is much shorter 
than any length scale of interest. For the protogalaxies that 
we study, this condition is generally satisfied; for example, 
the mean free path of a photon just above the Lyman limit 
in the core region of a z = 10 protogalaxy is only 0.01 pc 
(and is even smaller in higher redshift objects). It breaks 
down somewhat in the outermost regions of our protogalax- 
ies, where the density is several orders of magnitude lower, 
but H 2 formation in these regions is negligible. 

These approximations allow us to write the radiative 
transfer equation for an arbitrary line of sight through the 
protogalaxy 



dl v 
ds 



= -Kv{s)I v + jv(s), 



(38) 



(where I v is the specific intensity, k v and j v are the opacity 
and emissivity respectively, and where s is the position along 
the line of sight) in the simpler form 

dIv - _ ( \ T 
ds 

which can be further simplified to 

dlu 
dT„ 



-Iu. 



(39) 



(40) 



Solution of this equation is trivial; if the specific intensity 
at the edge of the protogalaxy is i„(0), then at a distance s 
into the protogalaxy, it is 



/„ = 7„(0)e- 



,0) 



(41) 



To obtain the mean specific intensity, we simply integrate 
this expression over solid angle. Using the conventional vari- 
ables r and /i (see figure |l|), we find that 



exp| 



(42) 
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Figure 1. A schematic diagram of the radiative transfer geometry 
used in our simulations. For each of the Ng mass shells, we pick 
a single central point (r^ and ry above are two examples) and 
compute column densities along N a lines of sight to this point. 
With each line of sight, we can associate an angle; for instance, 
with the j-th line of sight from )-&, we associate the angle Ojk- 
However, in most cases it proves simpler to work in terms of the 
direction cosine = cos0j k . 



where J v (r t ) is the incident intensity at the edge of the pro- 
togalaxy (i.e. the background radiation) and where t(v, r, fi) 
is the optical depth between r and the edge of the proto- 
galaxy in the direction of \i. 



3.4-1 Photoionization 

For the photoionization reactions (numbers 22-24, 26-27 & 
29), the opacity is dominated by continuum absorption by 
neutral hydrogen and helium. For static gas, we can write 
this as 



T{y,r,n) = a H (iy)N H (r,fi) + a U oN Hc (r,fj,) 
Y 



OH + 



4-4Y 



(43) 
(44) 



where the second line follows if we assume primordial abun- 
dances and a small fractional ionization (which also allows 
us to neglect He + absorption). 

Substituting this expression into equation W2l we obtain 



J<A r ) = 2^"( rt ) exp 
where 



Y 



4-4Y 



A/" H (r)= / exp[-JV H (r,/x)] d^i. 



(46) 



Aside from Afn(r), all of the terms in equation |45| are known 
at the beginning of the simulation. Moreover, if we assume 
that J v {rt) remains constant for the duration of the simula- 
tion, then it becomes possible to calculate the photoioniza- 
tion rates in advance, tabulating them as functions of 7Vh- If 
we do this, then instead of solving directly for the set of pho- 
toionization rates at each timestep, we can simply calculate 
Afn(r) and then interpolate the rates from the appropriate 
tables; this avoids a large number of frequency integrations, 
and results in a significant speed-up. 



3-4-2 Photodissociation 

For Lyman- Werner photons, the situation is a little more 
complicated, due to the large number of spectral lines in- 
volved. Fortunately, there is a simple approach that we can 
use which significantly reduces the amount of computation 
required to obtain the photodissociation rate. 

We begin by writing the total H2 photodissociation rate 
as a sum of the contributions from each individual rotational 
and vibrational level: 



fcdis«.H 2 



3 



k 7 Tb 7° 



(47) 



where hh 2 is the total number density of H2 molecules, 
and rij is the number density of molecules in the j-th ro- 
vibrational level. For each level, we can write the photodis- 



sociation rate in the form (Draine & Bertoldi 1996) 



^ J : ' " <1.V., . = 



dfi, 



(48) 



where Wij is the dimensionless equivalent width of the i- 
th Lyman- Werner line originating from level j, A r H 2 ,j is the 
column density of H2 molecules in level j, and /di s ,i is the 
dissociation fraction for level i (i.e. the fraction of excitations 
that are followed by dissociation, rather than de-excitation 
to some other bound state). 

The equivalent width can be written as 



[1 - exp(-a i (i')iVH 2 ,j)] — , 



(49) 



where the radiative cross-section cn(y) can be written in 
terms of the line profile cj>(i>, Vi) and oscillator strength / OSCi j 
as 



fa, 



(50) 



In general, the form of the line profile will depend upon the 
details of the temperature and velocity variations along the 
line of sight. For static, isothermal gas, however, it reduces 
to the familiar Voigt profile. Consequently, (?i{v) is known a 
priori, allowing us to tabulate Wij, and hence AWij/ dAfH 2 ,j, 
as functions of A^H 2 ,j at the beginning of the simulation. 
Similarly, we can tabulate the photodissociation rates kj as 
a function of 



)]A/H(r), (45) A/H 2lj (r) = 



dWi- 



dNi 



H 2 ,j 



dfi, 



(51) 



since the dissociation fractions /dis.i are determined purely 
by atomic physics and are independant of the physical con- 
ditions in the gas. In this way, we reduce the problem of 
determining the photodissociation rates to the simpler one 
of determining the column density of H2 molecules in each 
ro- vibrational level for a large number of lines of sight. Given 
a set of column densities for a particular level, we can obtain 
the corresponding values of dWij / dNu 2 ,j by interpolation; 
from these, we can derive Afn 2 ,i( r ) by numerical integration; 
and finally, the photodissociation rate kj is obtained by an- 
other interpolation. 

Two potential problems remain with this method, how- 
ever. Firstly, equation ^9] breaks down at high column den- 
sities as two effects that we have neglected - absorption 
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by neutral hydrogen, and overlap of the individual Lyman- 
Werner lines - become significant. However, these effects are 
important only for H2 column densities Nn 2 > 10 20 cm -2 ; 
since the largest H2 column densities that we encounter 
in our simulations are typically no larger than 7Vh 2 = 
10 18 cm -2 , it is generally safe to ignore these complications. ^| 

The second and more significant problem is the fact 
that, in the method as outlined, we have to deal with all of 
the rotational and vibrational levels of H2. This makes cal- 
culating the photodissociation rate much more costly than 
calculating the photoionization rates, as we need to com- 
pute several hundred different column densities for each line 
of sight. Similarly, the chemical network is also greatly en- 
larged and consequently takes longer to solve. Moreover, 
the excited vibrational levels typically have characteristic 
timescales of the order of 10 6 s, much smaller than the 
timescale on which the total H2 abundance changes, forc- 
ing us to solve the chemical network with a timestep that is 
orders of magnitude smaller than we could otherwise use. 

To overcome these problems, we make another approx- 
imation. Since, at the temperatures and densities encoun- 
tered in our simulations, most of the excited levels have very 
small populations, their net contribution to the photodisso- 
ciation rate is also very small. We can therefore safely ignore 
them, and concentrate only on the few rotational levels that 
are responsible for the bulk of the photodissociation. 

The simplest approximation that we can make is to as- 
sume that all of the H2 molecules are in the form of para- 
hydrogen and are to be found i n the ground stat e. This ap- 
proximation was suggested by Abel et al. (1997), and sup- 
poses that any ortho-hydrogen formed will be rapi dly con- 



verte d to para- hydrogen by collisions with protons (Gerlicb 



199C) 



H 2 (ortho) + H+ -> H 2 (para) + H H 



(52) 



Unfortunately, this is only valid at very low temperatures, 
typically T<50K. At higher temperatures, conversion from 
para-hydrogen back to ortho-hydrogen via the inverse re- 
action becomes competitive. Moreover, at realistic temper- 
atures and densities, ortho-para interconversion via proton 
collisions with excited states will also occur, and in general 
the ortho to para ratio will be close to its equilibrium value 
of three to one. 

We can therefore obtain a much better approxima- 
tion to the true photodissociation rate by including ortho- 
hydrogen. For simplicity, we assume that the ortho to para 
ratio is fixed and is exactly three to one. We continue to 
ignore the other excited rotational levels. 

Comparing the results of simulations performed using 
this approximation with those of simulations in which the 
rotational and vibrational level populations are set to their 
local thermodynamic equilibrium (LTE) values, we find that 
it performs well: although differences in the final H2 abun- 
dances are often observed, the effect on T CT i t remains small 
- it is increased by at most 30-40% at any given redshift, 
and by much less than this in the limits of low and high 
redshift. The uncertainty that this introduces is comparable 
to that due to our poor knowledge of the high-redshift ultra- 



violet background (see figure ^| below). Moreover, since the 
true level populations will lie somewhere between our simple 
approximation and the LTE case, this suggests that our ap- 
proximation is sufficient for our current purposes, and that a 
more detailed (and hence slower) treatment is not required. 



3.4-3 Calculating the column densities 

As demonstrated in the previous sections, we can determine 
the photochemical reaction rates for a given mass shell by 
computing the functions 



exp[-JV H (r,/Lt)] d/z, 



L ?*• 



AW, 



dNi 



— dfx, 



H 2 ,j 



(53) 



(54) 



A/k(r) 
and 

■A/k 2>j (r) = 

(where j = 0, 1 for the para and ortho-hydrogen ground 
states respectively) and then interpolating the rates from 
a pre-built look-up table. Moreover, knowledge of Nn(r,/j,) 
and Nn 2 ,j(r,n) is sufficient to specify Nb.(v) and A/k 2 ,j('") 
respectively, so the problem is essentially reduced to one of 
calculating these column densities. 

To do th is , we use a method very similar to that in 
Kepner et al. (1997). For each mass shell k, we choose a 
point Tk at the centre of the shell (see figure [l]). For each 
of these points, we compute the column densities along N a 
lines of sight with uniform angular separation. The column 
density of a chemical species i along one of these lines of 
sight is given by 



Ni(r,fj,)= / m[r'(x)]dx, 



(55) 



where L = rfi+[r 2 — r 2 (l — fi 2 )} 1 ' 2 is the total distance along 
the ray to the edge of the protogalaxy, and where r'(x) = 
(r 2 + x 2 — 2rxfi) 1 ^ 2 is the distance from the point labeled by 
x to the centre of the protogalaxy. Solution of this equation 
by means of numerical integration is straightforward. 

Given the set of column densities for a particular point 
rfc, we use equations [53] and ^ to compute A/h (?**;) and 
A/H 2 j(rfc), ^ and then finally use these values to compute 
the photochemical rates at that point. Finally, the spherical 
symmetry of the problem allows us to generalize these rates 
to the whole of the mass shell containing rjfe. 

The accuracy of this procedure clearly depends both 
upon the spatial resolution of the grid (N g ) and upon the an- 
gular resolution with which we sample it (N a ). However, as 
the density, temperature and chemical abundances all vary 
smoothly within the protogalaxy, we find that very high res- 
olution is not required: setting N g = 100 and N a = 20 proves 
to be sufficient for accurate calculation of the photochemi- 
cal rates. Simulations run with substantially higher values 
produce very similar results. 



3 In any case, it is possible to reformulate equ a tion | 49 so as 
avoid these difficulties, as in Draine & Bertold: (|1996|). 



to 



Tn the latt 



et al. 



(L993a, 



tions from Abgrall et al. (EOO0) 



Abgral] 



case we make use of molecular data from 
and iRoucfil ( 1997 ) , together with dissociation frac- 
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Process 



Reference 



Atomic line cooling: 
Hi 
He i 
Hen 

Molecular line cooling: 
H 2 

Collisional ionization: 
Hi 
He i 
Hen 

Recombination: 
Hn 
Hen 
He III 

Other chemistry: 
H — formation 



Bremsstrahlung: 

Compton scattering: 

Photoionization: 
Hi 
He i 
Hen 

Photodissociation: 
H 2 




|Brav et al\ (]2000fl 
Aggarwal et al (|L992|) 



Le Bourlot et al. (1999) 



Janev et al 



Jancv et al 



?\bcl et al. 



(1987) 



(1987) 



(1997) 



Fcrland et al. 



(1992) 



Aldrqvandi fc Pcauignoq ( 1973 ) 
Ferland et al.lnm.' 



Shapiro fe Kang ( 1987 ) 
Shapiro & Kang (1987) 



Spitzer (1978) 



Peebles 



(1993) 




(199J) 



Black fe Dalgarnc (1977) 



Table 2. The main processes responsible for heating and cool- 
ing metal-free gas, together with references to the source(s) from 
which the heating and cooling rates used in our simulations were 
taken. 



3.5 Heating and cooling 

The heating and cooling processes included in our model 
are listed in table ^, together with references to the sources 
for the rates adopted. At low temperatures, cooling is dom- 
inated by rotational and vibrational line emission from Eb; 
we include th is acco rding to the recent prescription of [La 



Bourlot et al. ( 1999 ), which is probably the most accurate 
available to date. We do not include the effects of HD cool- 
ing: this becomes significant only at very low temperatures, 
and in any case does no t appear to have a great effect on 
the evolution of the gas (Bromm et al. 2002). 

At high temperatures (T > 8000 K), cooling through 
the atomic lines of hydrogen (often referred to simply as 
Lyman-a cooling) rapidly becomes dominant. This sets an 
upper limit on T cr i t ; for simplicity, we neglect any redshift 
dependence of this limit, and assume that protogalaxies with 
T v i r ^ 8000 K can always cool successfully. 

Heating of the protogalactic gas is driven primarily by 
the photoionization of hydrogen and helium. H2 photodisso- 
ciation also contributes to the heating rate, but is generally 
not significant, due to the small H2 abundance. 



3.6 Running the simulation 

Given the initial temperature and chemical abundances, plus 
the set of chemical reaction rates, actually solving for the 
thermal and chemical evolution is relatively easy. As in the 
IGM case, we simply solve the coupled set of chemical rate 
equations 

dn 



— rr = Ci(rij,T) - Di(nj,T)m, 



dt 

AT _ 2 
dt 3fcn to t 



(r — A) , 



(56) 



(57) 



using the STIFBS integrator of Press et al. (1986) 



At the start of each timestep (hereafter time i), we com- 



pute the photochemical rates as outlined in sections 3.4 



We then use STIFBS to solve for the new chemical abun- 
dances and new temperature at the end of the timestep (time 
t + At), repeating this for each shell in turn. We next store 
these values, return to time t, and recalculate them using 
the same procedure, but with two timesteps of length At/2. 
We recalculate the photochemical rates at the intermediate 
point. 

We next test for convergence by comparing our two sets 
of results. If any of the chemical abundances or temperatures 
of any of the shells differ by more than 0.1%, then we re- 
ject the results and begin again from time t with a smaller 
timestep. Otherwise, we check to see whether we need to 
halt the simulation, using the criteria discussed below, and, 
if we do not, we begin the computations for a new timestep 
starting from t + At. 

One final approximation that we find useful in practice 
is to fix the H~ and H^" abundances at their equilibrium val- 
ues. This allows the integrator to take much larger timesteps 
than would otherwise be possible, but introduces very little 
error into the computed H2 abundances. 



3.7 Halting the simulation 

Using the method outlined in the preceding sections, we 
compute the chemical and thermal evolution of the proto- 
galaxy until one of two conditions is met: either the proto- 
galactic gas begins to cool strongly, or we exceed a preset 
time limit, tii m . 

To assess whether gas cooling is 'strong' enough requires 
an objective cooling criterion. A number of different possibil- 
ities have been sugge sted in the literatur e ( | Rees fc Ostriker 



1977; Iegmark et al. 



1997; Haiman et al 



20001). In. 



lations, we adopt the criterion used by |Haiman et al. (2000) 



we require that the elapsed time exceeds the cooling time, as 
computed at the edge of the protogalactic core, at a distance 
ro = 0.034r t from the centre of the protogalaxy. The advan- 
tage of this choice is that it avoids giving us a false positive 
result in cases where i coo i drops briefly below td yn at early 
times, but remains so for a time t <C t C ooi- As an additional 
sanity check, we also require that the final temperature be 
smaller than the initial temperature. 

If the protogalactic gas does not cool strongly, then the 
simulation will terminate when it reaches tn m . This pre-set 
time limit is required on purely practical grounds, to prevent 
simulations in which the gas does not cool from running for 
excessive amounts of time, but also has a physical justifica- 
tion. In our simulations, we treat protogalaxies as isolated 
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objects, uninfluenced by external events. In reality, they are 
part of a dynamically evolving mass distribution, and the 
majority will only survive for a limited time before merging 
with other protogalaxies of a similar or larger size. It is pos- 
sible to use the Press-Schechter formalism to calculate the 



distribu tion o f survival times as a function of mass (Lacey 
& Cole| [l993 ) , but for rare objects the mean survival time 



is typically of the order of the Hubble time and thus for 
simplicity we set ti im = tn- 



4 RESULTS 

In the following sections, we present results from a number 
of simulations that examine the effects of the X-ray back- 
grounds prod uced by th e various source models discussed in 
sections 



2.1.3 and 2.1.4 



Unless otherwise noted, all of these 
simulations assume the same parameters for the star forma- 
tion model: a standard Salpeter IMF, with M m i n = 0.1 Mq 
and A/ max = 100 Mq, a star formation efficiency e — 0.1 and 
a starburst lifetime t on — 10 7 yr. Additionally, all of the sim- 
ulations use the s ame cosmological model, the ACDM con- 
cordance model of |Wang et al. ( 200C ) , which has parameters 
(fi A , «m, ^6, h, n, cr 8 ) = (0.67, 0.33, 0.041, 0.65, 1.0, 0.9). 



4.1 Simulations without an X-ray background 

Before examining the effect of an X-ray background on pro- 
togalactic evolution, we first briefly study the evolution of 
T cr it in its absence. As well as allowing us to determine 
the sensitivity of our results to variations in the UV source 
model, this also provides us with a necessary baseline against 
which to compare our other results. 

In figure ^, we plot the evolution of T cr it with red- 
shift for t he thr ee different star formation models discussed 
in section 2.1.2. All three models assume a Salpeter IMF, 

= 100 M , as well as 



with maximum stellar mass M max 
our standard star formation efficiency, starburst lifetime and 
cosmological model, described previously. Our basic model 
assumes a metal-free stellar population, with a minimum 
stellar mass of M m i n = 0.1 Mq; the corresponding results 
are given by the dotted line in figure ^[ The dashed line 
illustrates the effect of increasing the minimum mass to 



M n 



1 Mq; the solid line assumes the same M n 



to- 



gether with a metallicity Z — 0.02Zq. 

Figure |^ demonstrates that although the strength of 
the Lyman- Werner background increases by almost a fac- 
tor of five as we move from our basic model to the metal- 
enriched model, this has less effect on T cr i t than we might 
expect: the difference in T cr i t between the three models is 
never more than 50%, and the qualitative details of its evo- 
lution are very similar in all three models. This suggests 
that the uncertainty introduced by our lack of knowledge 
of the properties of the primordial stellar population need 
not be unduly limiting. However, it is also clear that some 
uncertainty remains, and this will place a lower limit on the 
magnitude of any effect that we can reliably claim to detect, 
as small variations in T cr it due to the X-ray background will 
be swamped the by the error resulting from the uncertainty 
in the Lyman- Werner background. 

In the work that follows, we take as our baseline the 
results of our basic, metal-free, low M m i n model; as figure ^ 




Figure 2. The evolution of the critical temperature as a func- 
tion of redshift in the absence of an X-ray background. Results 
are plotted for three different UV source models: two metal-free 
models, one with a minimum stellar mass M m ; n = 0.1 Mq and the 
other with M m ; n = 1 Mq (represented by the dotted and dashed 
lines respectively) and a Z = 0.02Zq model with M mln = 1 Mq 
(solid line). In all three cases, the IMF is of Salpeter form, with 
a maximum mass M max = 100 Mq. The remaining parameters 
of the star formation model (star formation efficiency and star- 
burst lifetime) are the same in all three simulations, as is the 
cosmological model. 



demonstrates, this minimizes the strength of the Lyman- 
Werner background, and thus will tend to maximize the ef- 
fectiveness of the X-ray background. 



4.2 Massive X-ray binaries 

In figure ^, we plot the evolution of Tcrit in the presence 
of the X-ray background generate d by the massive X-ray 
binary model described in section 2.1.3. For the purposes 
of comparison, we also plot the results of our basic X-ray 
free model. Initially, the evolution of T cr i t is the same in 
both models, implying that the X-ray background has little 
or no effect on the gas. At a redshift z ~ 20 and critical 
temperature T cr it — 2000 K, however, the models begin to 
diverge significantly. In the X-ray free model, the critical 
temperature continues to increase rapidly until it reaches 
its maximum value of T cr i t = 8000 K. In the X-ray binary 
model, by contrast, the rate of increase of T cr i t is significantly 
slower, and it fails to reach its maximum value by the end 
of the simulation at z — 10 Ln 

It is clear from figure plthat the presence of the X-ray 
background significantly affects the evolution of gas in the 
larger of the H2-cooled protogalaxies. In small protogalax- 
ies, on the other hand, the negative feedback caused by the 
Lyman- Werner background remains as strong as ever. 

A simple way to judge the importance of this effect is 
to examine the difference it makes to the fraction of gas 



5 We choose to end our simulations at z = 10 because at lower 
redshift we expect the effects of ionizing radiation from both stel- 
lar sources and quasars to become increasingly important and 
thus our results to become unreliable. Clearly if reionization oc- 
curs at z > 10, the same is true for some portion of the results 
plotted here. 
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Figure 3. The evolution of T cr i t as a function of redshift for an 
X-ray bac kgrou nd produced by the massive X-ray binary model 
of section 2.1.3 (dashed line). For comparison, we also plot the 
evolution of T cr ; t for our basic X-ray free model (solid line). 



Figure 4. The ratio of the cooled gas fraction in the massive 
X-ray binary model to that in the X-ray free model. The mass 
fractions are calculated using the Press-Schechter formalism, as 
outlined in the text. 



in the unive rse that can collapse and cool. As we saw in 
we can use the Press-Schechter formalism to 



section 



2.1.1 



express the cool gas fraction as 



F(z, M crit ) = erfc 



Sc(z) 



V2a{M clit ) 



(58) 



Using the relationship between M cr it and T cr i t derived in 
that section, it is straightforward to calculate the evolution 
of F(z, M cr it) in both the X-ray binary and X-ray free mod- 
els. 

To better highlight the difference between the two mod- 
els, we plot in figure [l| the ratio of the cooled gas fraction in 
the X-ray binary model to that in the X-ray free model. At 
high redshift, the evolution of T cr it is the same in both mod- 
els, and thus the ratio is one. Below z = 30, the behaviour of 
the models begins to diverge, but this has little effect on the 
ratio until the models begin to diverge sharply at z ~ 20. 
Thereafter, the ratio rises sharply to a peak at z ~ 15, and 
subsequently declines as the growth in the cooled mass frac- 
tion becomes dominated by the formation of larger proto- 
galaxies that cool via Lyman-a emission. 

Figure |^ shows us that by ignoring the effect of the 
X-ray background we underestimate F(z, Merit) by at most 
a factor of two. In fact, the difference is likely to be even 
smaller: we have assumed that all of the gas in a proto- 
galaxy with M > M cr it will cool, but this need not be the 
case if an X-ray background is present, as X-ray heating will 
prevent gas from cooling in the low-density outer layers of 
the protogalaxy. Whether this will affect the star formation 
rate within individual protogalaxies is not clear, but in any 
case the difference in the globally-averaged star formation 
rate is unlikely to be greater than a factor of two. Note, 
however, that this extra star formation occurs entirely in 
low-mass systems, from which ion izing photons (|Ricotti & 
Shul 2000) and supernova ejecta (Ferrara & Tolstoy 200C) 
can more readily escape. Doubling the star formation rate 
may thus significantly increase the feedback of primordial 
star formation on the IGM. 

Finally, given the uncertainties in the data und erlyin g 
our simple massive X-ray binary model (see section 2.1.3), 




Figure 5. As figure ^, but including variants of the X-ray binary 
model with 10% (dash-dotted line) and 1000% (dotted line) of 
the flux of the basic model (dashed line). As before, we also plot 
the evolution of T cr it for the X-ray free model (solid line). 



to changes in the strength of the X-ray background. Ac- 
cordingly, we have run additional simulations in which the 
strength of the X-ray sources was increased or decreased by 
a factor of ten. The results are plotted in figure ^, together 
with the results of our basic X-ray binary model and of the 
X-ray free model. 

Unsurprisingly, increasing or decreasing the strength of 
the X-ray background alters the evolution of Tcrit. Decreas- 
ing it by an order of magnitude increases T cr it to the point 
where its evolution is little different from that in the X-ray 
free model. Increasing it by an order of magnitude, on the 
other hand, systematically lowers Tcrit, although never by 
more than a factor of two. 



4.3 Supernova remnants 



As we discussed in section 2.1.4, it is possible that mas- 



it is of interest to investigate the sensitivity of our results 



sive X-ray binaries are much less abundant at high redshift 
than at the present day. If so, then the X-ray emission of 
star-forming galaxies will be dominated by supernova rem- 
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Figure 6. As figure ^, but for an X-ray background produced 
by inverse Compton emission from supernova remnants. We as- 
sume that 10% of the total supernova energy is transferred to the 
relativistic electrons powering the emission, and again plot the 
results of the X-ray free model for comparison. 



O 
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Figure 7. The evolution with redshift of the fractional abun- 
dance of H2 in the intergalactic medium, plotted for the X-ray free 
model (solid line) and the massive X-ray binary model (dashed 
line). 



In this con text, the possible observational tests suggested by 



Oh 



(2001) could prove extremely valuable. 



nants (SNR). These will generate X-rays through two main 
emission mechanisms: thermal bremsstrahlung from hot gas, 
and inverse Compton scattering of the CMB by relativistic 
electrons. We consider the effects of these mechanisms sep- 
arately. 

We examined two possible variants of the thermal 
bremsstrahlung model. In one, we assumed that the char- 
acteristics of the emission are broadly the same as those 
observed locally, with a fraction / x = 2 x 10 -4 of the su- 
pernova energy being radiated as X-rays with a charac ter- 
istic temperature T K = 1 keV (Helfand & Moran 2001). In 
the other model, we assumed that all high-redshift super- 
novae explode in extremely dense surroundings, producing 
X-ray bright, ultra-compact remn ants with / x = 0.01 and 
T x = 30keV (fTerlevich et oDpgl). Realistic models should 



lie somewhere between these two extremes. 

However, we found that in neither of these cases does 
the X-ray background affect the evolution of T C rit: at our 
level of temperature resolution, the results are identical to 
those obtained for the X-ray free model. The obvious con- 
clusion is that the background produced by bremsstrahlung 
is simply too weak to be effective. 

Our inverse Compton model fares somewhat better. In 
this model, the strength of the X-ray background is propor- 
tional to the mean fraction f e of the supernova explosion 
energy that is transferred to relativistic electrons within the 
remnant. This value is not known accurately and so we con- 
sidered two possible cases, one with / c = 10 -3 and another 
with / e = 0.1, these being conservative lower and upper 
bounds on the true value. In the former case, we again saw 
no significant effect on the evolution of T cr i t . In the latter 
case, on the other hand, we saw results very similar to those 
obtained for the X-ray binary model, as illustrated in fig- 
ure u 

In principle, therefore, inverse Compton emission from 
supernova remnants could be as important an X-ray source 
as emission from X-ray binaries. Ultimately, however, its 
importance depends upon the value of / e , and no firm con- 
clusions are possible until this value is better constrained. 



4.4 The thermal and chemical evolution of the 
IGM 

As well as determining the evolution of T cr i t , our simulations 
also allow us to study the thermal and chemical evolutio n of 



the diffuse intergalactic medium, as outlined in section 3.2 
As a simple example, we plot in figure [?] the evolution of 
the fractional H2 abundance in the IGM for the X-ray free 
model (solid line) and the X-ray binary model (dashed line). 

In both cases, the fractional abundance rapidly de- 
creases from its primordial value due to photodissociation by 
the ultraviolet background, reaching /h 2 = 10~ 9 by 2 ~ 35. 
Subsequently, its decline slows, in part because the rate of 
increase in the strength of the ultraviolet background also 
slows. Below z — 20, the behaviour of the two models di- 
verges. In the X-ray free model, /h 2 continues to decline un- 
til the end of the simulation. In the X-ray binary model, on 
the other hand, the increasing ionization of the IGM boosts 
the H2 formation rate to the point where it overtakes the 
photodissociation rate and the H2 abundance, after reach- 
ing a minimum at 2 ~ 17, begins to climb. Nevertheless, it 
remains extremely small at the end of the s imulation, read- 
ily justifying our assertion in section [2.2.1 that absorption 
by intergalactic H2 does not play a significant role in deter- 
mining the strength of the Lyman- Werner background. 

In figure ^, we plot the evolution of the fractional ion- 
ization of the IGM. In the X-ray free model, this remains ap- 
proximately constant over the lifetime of the simulation, as 
the recombination timescale is significantly longer than the 
Hubble time. In the X-ray binary model, on the other hand, 
photoionization by the growing X-ray background eventu- 
ally overcomes the very small recombination rate and drives 
the fractional ionization upwards, increasing it by just over 
an order of magnitude by the end of the simulation. Even so, 
it remains small, demonstrating that the X-ray background 
does not contribute significantly to cosmological reioniza- 
tion. 
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20 

Redshift 




Figure 8. The evolution with redshift of the fractional ioniza- 
tion of the intergalactic medium, plotted for the X-ray free model 
(solid line) and the massive X-ray binary model (dashed line). 



Figure 10. The evolution with redshift of the temperature of the 
intergalactic medium, plotted for the X-ray free model (solid line) 
and the massive X-ray binary model (dashed line). 



20 

Redshift 



Figure 9. The evolution with redshift of the fractional ionization 
of helium in the intergalactic medium, plotted for the massive X- 
ray binary model. In the X-ray free model, the He+ abundance 
remains negligible throughout the simulation. 



The effect on the fractional ionization of helium (the 
ratio of He + to He) is rather more striking. Th e post 



reco mbinat ion He + abundance is extremely small (standi 
et al. 199S ), and in the X-ray free model remains at this low 



level throughout the simulation. In the X-ray binary model, 
on the other hand, it increases dramatically over the course 
of the simulation, reaching / Hc + — 4x 10~ 3 by z = 10. This 
is illustrated in figure H. 

Finally, in figure lid we plot the evolution with redshift 
of the temperature of the IGM. In the X-ray free model, 
adiabatic cooling dominates the thermal evolution, and the 
temperature falls off approximately as (1 + z) 2 . In the X- 
ray binary model, on the other hand, photo-electric heating 
begins to heat the IGM strongly at z = 20, driving the 
temperature up to T = 100 K by the end of the simulation. 

Thus, although the X-ray background does not con- 
tribute significantly to the reionization of the IGM, it does 
produce substantial reheating prior to reionization. More- 
over, given the large mean free path of the X-ray photons, 
this reheating occurs almost uniformly throughout the IGM, 



rather than being localized to the vicinity of star-forming 
galaxies. 

One consequence of this reheating is that the formation 
of very small-scale structure will be suppressed, as the in- 
creased temperature of the IGM leads to an increased Jeans 
mass. This is unlikely to affect the global star formation 
rate, however, as star formation within these small struc- 
tures would in any case be strongly suppressed by the ul- 
traviolet background. Nevertheless, it will reduce the mean 
clumping factor of the IGM below the level that we would 
otherwise pred i ct, which may in turn speed up reionization 
flHaiman et"al||200l| ; [Barkana fc Loebj ^002[ ). 

Reheating also affects the visibility of t he IGM in the 



redshifted 21 cm line of neutral hydrogen. Madau et al. 



( 1997 ) show that scattered Lyman-a emission from high red- 
shift galaxies efficiently couples the spin temperature of the 
Hlhyperfine levels to the kinetic temperature of the gas. If 
the gas temperature is smaller than the CMB temperature, 
this results in 21 cm line absorption; if it is greater, then 
it results in emission. Absor ption is easier to detect than 
emission (scott & Reef 1990), but the heating produced by 



the X-ray background implies that absorption occurs only at 
z ^ 15, and that concentrating on detecting 21 cm emission 
may be the more viable strategy. 



5 CONCLUSIONS 

The results of the previous section allow us to assess the im- 
pact of the high redshift X-ray background that is produced 
by star-forming galaxies. If we assume that the X-ray emis- 
sion of these galaxies is similar to that observed locally, and 
that the same correlation between X-ray luminosity and star 
formation rate applies, then we find that the background 
produced is strong enough to partially offset the effects of 
UV photodissociation in large (T v i r > 1000 K), H2-cooled 
protogalaxies. 

However, local emission is dominated by massive X-ray 
binaries, which may not form in large numbers at high red- 
shift. Therefore, we have also explored the effect of an X-ray 
background produced by emission from supernova remnants. 
If this emission is dominated by inverse Compton scatter- 
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ing and if the fraction of the supernova energy transferred 
to the relativistic electrons powering this emission is large, 
then the resulting background has very similar effects to one 
produced by X-ray binaries. On the other hand, if the frac- 
tion of energy transferred is small, then the background has 
little or no effect. 

In addition to inverse Compton emission, we have also 
examined the effect of thermal bremsstrahlung emission 
from hot gas in the remnants, and find that even if all super- 
novae were to form X-ray bright, ultra-compact remnants, 
the resulting X-ray background would still be too weak to 
significantly affect protogalactic evolution. 

Finally, none of these models produces an X-ray back- 
ground that is strong enough to balance UV photodis- 
sociation in small protogalaxies, with virial temperatures 
T v i r < 1000 K. In these protogalaxies, negative feedback al- 
ways dominates. 

How significant are these results? One simple way to 
assess this is to study the evolution of the mass fraction 
of cooled gas, which represents the total amount of matter 
available to form stars. Comparing its evolution in the X-ray 
binary model to that in the absence of an X-ray background, 
we find that it is increased by approximately a factor of two. 
Given our star formation model, this corresponds to an in- 
crease in the global star formation rate by the same amount. 
However, this is small compared to the order of magnitude 
increase that would result if we were simply to ignore the 



effect of the UV background (see figure 7 of Haiman et al 
200C|). 



In reality, the difference between the two models may 
be greater than this because the additional star formation 
takes place entirely in low-mass systems, from which ionizing 
photons and supernova-produced metals can readily escape 



(Ricottifc Shull 



2000 



Ferrara fc Tolstoy 2000). Doubling the 



star formation rate may therefore have more than double the 
impact on the intergalactic medium. However, to properly 
assess the ultimate importance of this effect requires more 
detailed modelling, beyond the scope of this paper. 

Ultimately, understanding the history of star formation 
in the small protogalaxies studied in this paper remains im- 
portant even if they do not contribute to the reionization 
or enrichment of the IGM to any great degree. This is sim- 
ply because, in a hierarchical universe, these protogalaxies 
are the building blocks from which larger galaxies form and 
therefore set the initial conditions for later stages of galaxy 
formation. In particular, very little metal enrichment of the 
primordial gas is required in order to allow the CNO cycle 
to operate and population II (rather than population III) 
stars to form, and yet this can have a profound effect on 
the predicted spectral energy distribution of an early stellar 
population. 

Although the main purpose of our study was to exam- 
ine the effects of the X-ray background on the thermal and 
chemical evolution of gas within protogalaxies, our approach 
also allows us to examine the effects of the background on 
the diffuse IGM. Our main results are three-fold: 

(i) We confirm the rapid destruction of H2 in the inter- 

~ ( |200C| ), but also 



Haiman et al 



galactic medium noted by 

show that when an X-ray background is present the H2 
abundance does not continue to decline indefinitely, but 
eventually stabilizes and may even begin to increase. How- 



ever, it never becomes large enough to significantly affect 
the Lyman-Werner background. 

(ii) We show that although photoionization by the X-ray 
background significantly increases the fractional ionization 
of the IGM (and in particular the fractional ionization of 
helium) , the bulk of the gas remains mostly neutral, demon- 
strating that the contribution of the X-ray background to 
cosmological reionization is small. 

(iii) We find that the X-ray background will also heat the 
intergalactic gas, raising its temperature to T ~ 100 K by 
z — 10 (compared to T = 3.5 K in the X-ray free model). 
This will suppress the formation of structure on the smallest 
scales by increasing the Jeans mass. It is unlikely to affect 
the global star formation rate, since T C rit S> 100 K, but may 
speed up the process of reionization by reducing the mean 
clumping factor of the IGM. 

A number of other authors have studied the effects of 
radi ative feedback on th e f ormation of Hg-cooled protogalax- 



(ICiardi et al\ k000a|,|bl; E aiman et al. 
" p00 412002^ 



1996b 



1997, 



200C 



Machacek et al [2001 , |2002[ ). Most of these studies have con- 
centrated solely on the effects of the Lyman-Werner back- 
ground, generally finding that it suppresses cooling (and 
hence star formation) by z ~ 20-30, prior to cosmological 
reionization. Our X-ray free simulations support this conclu 



sion. I n par ticular, comparison with the results of Haiman 



et al 



(200C), who use a very similar method but with a dif- 
ferent implementation of radiative transfer, gas chemistry 
and cooling, shows good agreement. 

Rather less work has be en done on the effects of the high 
redshift X-ray background. Haiman et al. ( 1996b ) examined 
the formation of molecular hydrogen in a constant density 
primordial cloud illuminated by a power-law UV spectrum 
extending into the hard X-rays, and found that in some 
case s, the background c ould enhance H2 formation. This 



lead 



Haiman et al. 



(200C ) to examine the effects of a quasar- 
produced X-ray background, using a very simple model in 
which the X-ray flux is a fixed fraction of the Lyman-Werner 
band flux, modulated by absorption by a fixed column den- 
sity of H and He. They found that an X-ray to UV flux ra- 
tio of 10% was enough to overcome negative feedback, and 
that a higher flux r atio could potentially produce positive 
feedback. Recently, (Machacek et al. 2002) have studied the 
effects of similar model backgrounds using an adaptive-mesh 
hydrodynamics code. They also find that X-rays reduce the 
effectiveness of negative feedback, but that the latter still 
dominates. They do not find evidence for the positive feed- 
back predicted by Haiman et al. (200C). However, their sim- 



ulations do not include the effects of H2 self-shielding, and 
thus potential ly underestimate the a mount of H2 that forms 



(although see 



Machacek et al. 



2001 



for a different view). 



There are several significant differences between our 
work and these previous investigations. Firstly, we do not as- 
sume a fixed spectrum or intensity for the X-ray background; 
rather, we specify the properties of the X-ray sources and 
subsequently compute the build-up of the background in a 
self-consistent fashion. Moreover, we consider source models 
where the X-ray emission is proportional to the star for- 
mation rate, as is observed to be the case for star-forming 
galaxies at low redshift; the relationship between star for- 
mation rate and X-ray emission in quasar-based models is 
far less clear. As a result, we generally consider X-ray back- 
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grounds signficantly fainter than those studied in the papers 
cited above. This makes direct comparison of ou r result s dif- 



ficult. However, we note that Machacek et al. (2002) find 
X-ray feedback to be ineffective below M v i v = 10 Mq, re- 
gardless of the strength of the X-ray background; at the 
redshifts they consider, this corresponds to a virial temper- 
ature TUr ~ 1000 K, and thus agrees well with the similar 
result obtained in this paper. 

In closing, we note that a number of uncertainties still 
remain in our treatment of this problem. Some of these - the 
high redshift star formation rate or the appropriate popu- 
lation III initial mass function, for instance - we simply do 
not know at the present time. Fortunately, changes to our 
assumed values can be readily incorporated in the frame- 
work laid out in this paper. Other uncertainties arise from 
our method of simulation; in particular, from our assump- 
tion of a static density profile for the protogalactic gas. We 
hope to address these issues in future work. 
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